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ABSTRACT 


The elementary problem considered was: Given earth- fixed space 

directions from two non-intervisible ground stations A and B to two elevated 
points 1 and 2, find the direction AB. Since more than two elevated points 
were to be used in the practice, a preliminary adjustment for each side 
AB was carried out. As the result of this adjustment , the two correlated 
quantities S ^ and t^, which were obtained, can be used as fictive 
observations , in the main adjustment of a network of several stations. 

If there are N stations, the main adjustment contains 3N parameters. 

Of these , four cannot be determined by triangulation. It hardly seems 
worthwhile, however, to include in the main adjustment such observations 
as satellite distances, ground tri angulations, levelings, and astronomical 
coordinates (with gravimetric deflections of the vertical) . It is suggested , 
therefore, that four other conditions be used, in order to keep the gravity 
center of all stations and the mean scale of the approximate coordinates 
fixed. In this way, the overall accuracy obtainable by the adjustment of 
the network can be estimated. 




1. The Fundamental Problem 


We start from the well-known elementary problem of intersection on 
a plane (see Fig. 1): 

Given the coordinates x, y for the corner points Al, A2, Bl, and B2 
of a quadrilateral, 

Wanted the coordinates for the intersection point Q of the diagonals 
Al - Bl and A2 - B2. 

The solution can be obtained by the following sequence of formulas: 


K 1 = X A1 y Bl - *B1 y . 


Al 


< L 1 = y Al ~ y Bl 

M 1 13 *B1 ' X A1 

/ 


K 2 = X A2 y B2 ' ^2 y A2 
'( L 2 = y A2 y B2 


**2 = ^2 " X 


A2 


N * L 1 M 2 - L 2 M 1 


/ x Q = (M 1 K 2 - M 2 K^/N 


y Q = (K 1 L 2 " K 2 L l )/N 


( 1 ) 


( 2 ) 


Suppose now that A and B are two observation stations between which 
there is a forest or a hill which makes a direct sighting from A to B 
impossible. Furthermore, suppose that, on the top of the obstacle, there 
are two fixed target points, 1 and 2, which can be identified and sighted 
both from A and from B. 

Let a and b denote the horizontal and vertical bearings, respectively, 
observed at the ground stations toward the elevated targets. Assume that 
the distance AB is so short that the curvature of the earth can be 
neglected and that the horizontal bearings, observed by a compass, are 
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Fig. 1. Intersection of diagonals of a quadrilateral on the 
celestial sphere gives the earth-fixed space direc- 
tion between two stations A and B at which two 
satellite positions 1 and 2 have been observed. 



referred to the same vertical zero plane* It is advisable to select a 
zero plane which is close to direction AB. 

Assume that, in the infinity behind station B, there is a vertical 
screen perpendicular to this zero plane. Then the bearings from A to 1 
and 2 give two points A1 and A2 on the screen. The bearings from B to 1 
and 2 are taken in reverse, replacing a by 180° + a and b by -b. Then two 
more points, Bl and B2, are obtained on the screen. The geometric solution 
of the fundamental problem now gives point Q, or the bearing from A to B. 

For the numerical computations, we must consider the screen as a 
gnomonic projection of the celestial sphere to the tangent plane at zero 
point: 


x = tan a 

(3) 

t an b 

y =5 

cos a 


Substituting these coordinates into (1) , 
reverse solution of (3) gives and b^. 


formulas (2) give x^, y^ and the 


2. The Equatorial Coordinates 

In the satellite traingulation, the horizontal coordinates a and 
b are replaced by equatorial coordinates a and 6. However, the zero plane 
of a is revolving with respect to the ground stations. Therefore, we 
introduce the Greenwich hour angles 


t » 0 - 


a 


( 4 ) 


where 8 is the Greenwich sidereal time of the observation. 

Because the satellites are moving targets, the observations at both 
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ground stations mist be made simultaneously. This can be accomplished by 
photographing flashing satellites against the background of stars — or 
chopping the trail of the satellite by synchronized shutters at both 
stations . 

Without going into the technical details of the measurements and their 
reductions , we suppose here that the positions t^ and 6^ of each event 
i at each station A, corrected for the refraction and other disturbing 
effects, are already available. In the same way the mean errors m fc and m^ 
should be readily estimated. 

Writing t instead of a and 6 instead of b, we may use formulas (3) 
and (1) - (2) for the calculation of t^ and 6^, which indicate the earth- 
fixed space direction from A to B. However, the following sequence of 
formulas is more convenient for the practical computations: 

For i = 1 and i = 2 , compute 


<i * 

tan 

S Bi sitt l Ai 

- tan 

6.. sin 
Ai 

fc Bi 


*i- 

tan 

S A1 cos Hi 

- t an 

S B1 COS 

fc Ai 

(5) 

p i = 

sin 

(t Bl - C Ai^ 






V 


Our numerical example has been taken from VerOffentlichung des Geod. 
Inst. Potsdam, Nr. 29, (1965) referring observations of Echo I in Potsdam and 
Bucharest. The reduction to the center of the balloon has been Computed 
again for each observation, instead of the mean of each series, and the 
results differ slightly from those published by Arnold and Schoeps (1965). 
Similarly, the weights of observations for our adjustment do not agree 
with those used in Potsdam. 


Observations: 

t 

6 

tan 6 


A1 

-27° 16 M4V70 
+13° 2 7* 18V 76 
0.23925217 


B1 

+ 5° 35 1 53V 50 
+34° 3.1' 11'.'45 
0.68779110 


A2 

-58°40'00V28 

+24°51 , 40'.'66 

0.46336373 


B2 

-39 o 05'57V90 
+57°16 '22V22 
1.55603704 
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sin t 
cos t 


- .45819587 
+ .88885125 


+ .09755152 
+ .99523048 


- .85415716 
+ .52001495 


- .63066790 
+ .77605283 


i 

K 

X 


cot 


sin 

tan 


M 

t. 


1 

-0.33848245 

-0.37323293 

+0.54271928 

-0.26479004 

-75°10'08V45 

-0.96668510 


2 

-1.03687155 

-0.44956779 

+0.33491434 


1x2 

+0.11898805 « v 

-0.44936755 

-0.23482380 


-0.50515590 
-26 8 48 '03721 


For the adjustment (see pages 14 and 15) , these approximate values 
have been rounded off to 


t = -75°10'10" and 6. - -26°48'00" 

Q Q 

In Arnold and Schoeps (1965) , there are 18 satellite positions (numbered 
1 - 18) which have been observed simultaneously at Potsdam and at Bucharest 
against the background of stars. From these double observations, numbers 1 
and 12 have been used for the approximate computation before the adjustment 
of all 18 events. 

3. The Preliminary Adjustments 

In the practice, more than two events are observed for the determination 
of the direction of each side AB. These observations usually are independent 
of those of other sides, e.g. , AC or BC. We omit here the slight correlation 
between sides AB, AC, and BC which occurs if some events are observed 
simultaneously at three stations A, B, C. 

For each side, we perform a preliminary adjustment. Before this 
adjustment, we need the first approximate values t^, 6^ which can be 
computed from two selected events as explained above or from approximate 
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For the estimation of the weights p^, we suppose that the mean error of 
one measured-event position is independent of the direction: 


”Ai ‘ m(S Ai ) * m(t Ai ) ’ c ° 8 6 A1 


but may be variable from one photograph to another. The variance of 1^ is 
then 


m. 


2 , 2 
m A , sin o A , 
Ai Ai 


t 2 . 2 
+ m^ sin a 


Bi 


sin^a 


AB 


(id 


where o means the arc from Q to Ai or Bi or from Ai to Bi, respectively 
The computation formulas read 


cos 0 M = sin 6^ sin 6^ + cos 6^ cos 6^ cos (t^ - t^) 


Si - 2 °AB - C ° s 2{ A1 C ° s2{ Bi (K i 2 + X l + "i 2) 


( 12 ) 


The weights p^ are, of course, inversely proportional to nu 


Figure 2 gives a graphical table of weights p in the case where m 


Ai 


m Bi' 


The normal equations of the adjustment read 




[paa] g 4- [pab] n + [pal] = 0 

[pab] £ + [pbb] n + [pbl] = 0 


( 13 ) 


After the solution of these equations, we have new approximate values 
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Safe! life 


9 



Fig. 2. Weight p of a satellite event as function of angles a = BAS and 




coordinates of stations A and B, as will be explained in the next section. 
However, the formulas (5) can be applied to all observed events i = 1, 2, 

3 ... n. Then formulas 

^i ^i 

tan t^ = — , tan 6^ = — cos t^ (7) 

K i K i 

give the coordinates t^, 6^ of that point on the celestial sphere which 
is at the pole of great circle Ai - Bi, The distance of point t^, 6^ 
from this great circle can be computed by formula 


1^ = p[sin 6^ sin 6^ 4- cos 6^ cos 6^ cos (t^ - t^) ] 


( 8 ) 


where p « 206265 seconds of arc. 

Differentiation of (8) shows that if 6^ and t^ are corrected by 
d6^ = g and dt^ = n/cos 6^, respectively, the new distances will be 


v. =a. £4b. n+1. 

l l l l 


(9) 


2 2 

where .a^ + b ” 1, approximately. 


In fact, we have 


a. - cos 6^ sin 6 , - sin 6^ cos 6, cos (t, - t rt ) 

x Q i Q i l Q 


b_^ = cos sin (t^ - t^) 


( 10 ) 


The least-squares solution is based on the condition 


[pw] = £ p^ v^ v^ - minimum 
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( 14 ) 



v 


+ £ 


tp + n/cos 6^ 


which are entered as fictitious observations into the main adjustment of 

a network. In addition, we need the weight matrix 

. \ 


1 


P 


AB 



[paa] 

[pab] 


V 

where m = 1 in ideal cases. 


[pab] 
[pbb ] 




) 

Altemat ive ly , 


(15) 


2 fpw] 

m = 

n - 2 


(16) 


can be determined separately for each side AB. Here n is the number of 
^events i used in the preliminary adjustment. 

The formulas (5) and (7) to (12) give the true geometric quantities 
if we want to illustrate the events graphically by the diagonals mentioned 
in the first section (see Fig. 3) . Especially 


a - cos 3 

( 

b - -sin 3 

where 3 is the angle between the diagonal and the meridian of point Q. 

For the programming of the numerical computations, however, the 
following sequence of formulas is more convenient: 
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- 75° 10 20 


- 75*1010 


- 75 ° 10 ' 0 " 


Error ellipse of the direction Pots dam-B ucharest, as 
computed from 18 satellite events. 




(Ql) 

* sin 

6 

Q 


(Q2) 

= cos 

6 

Q 



(Al) 

- sin 

6 ai 


(Bl) 

= sin 

<5„. 

Bi 



(A2) 

- cos 

6 a1 cos 

(t Ai t Q ) 

(B2) 

- cos 

6 B1 cos 

^Bi 

- V 

(A3) 

- cos 

6 ai sin 

(t A± " t Q ) 

(B3) 

" cos 

6 bi sin 


-v 

(A4) 

= (Ql) 

(Al) + 

(Q2) (A2) 

(B4) 

• (Ql) 

(Bl) + 

(Q2) 

(B2) 


K - (A3) (Bl) - (Al) (B3) 

M = (A2) (B3) - (A3) (B2) 

P = /mjj U = (A4) 2 } + n^ 2 {1 - (B4) 2 } 

L = 206265 {(Q1)M+ (Q2)K }/P 

A = { (Q2)M - (Ql)K } /P 

B = { (Al) (B2) - (A 2) (Bl) } /P 

2 

Normal equations: [A ] £ + [AB] n + [AL] - 0 

[AB ] K + [B 2 ] H + [BL] = 0 

[pvv] = [L 2 ] + ?[AL] + n[BL] 

The true distances of point Q from arcs Ai - Bi are 

L 

/a 2 + B 2 


Hence : 


A£ + Bn + L 

6F77- 
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Example 

> 3 

H* 

II 

J 3 

H- 

II 

H 4 




fc Q 

-75° 10' 10 '.'00 


°Q 

-26°48 , 00V00 

l Ai 

-2 7° 16 M4V70 


°Ai 

+13°27'18 , .'76 

Si 

+ 5°35'53V50 


°Bi 

+34°31 , 11 , .'45 

Si - 

t Q +47°53 ’55V30 

t 

Bi " t Q 

+80° 46 '03 ’.'50 

(Ql) 

-0.45087754 


(Q2) 

+0.89258582 

(Al) 

+0.23268520 


(Bl) 

+0.56669168 

(A2) 

+0.65204128 


(B2) 

+0.13219026 

(A3) 

+0.72159530 


(B3) 

+0.81325659 

(A4) 

+0.47709027 


(B4) 

+0.13751740 

K 

+0.21968928 


A 

+0.3679 

M 

+0.43488900 


B 

-0.2558 

P 

1.3242 


L 

+ 1 ”53 

Adjustment 




Operation equations 




i 

P A 

B 

L 

V 

1 

1.324 +0.3679 

-0.2558 

+1 "5 3 

-0V20 

2 

1.324 3668 

- 2582 

+0.41 

+2.31 

3 

1.323 3656 

- 2606 

+7.96 

-14.52 

4 

1.321 3626 

- 2666 

+1.25 

+0.39 

5 

1. 320 3613 

- 2690 

+4.95 

1 

00 

6 

1.206 +0.3459 

-0.2506 

-0.78 

+5.03 

7 

1.204 3436 

- 2544 

-2 . 14 

+8.19 

8 

1.203 3427 

- 2556 

-2.88 

+9.89 

9 

1.203 3419 

- 2569 

+1.37 

-0.04 

10 

1.139 +0.1959 

-0.3694 

+1 . 80 

-2.29 

11 

1.138 1930 

- 3704 

+1.03 

-0.48 

12 

1.287 +0.1431 

-0.4260 

+1.05 

-0.84 

13 

1.286 1406 

- 4259 

-1.81 

+5.50 
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14 

1.283 v 

1356 

“ 

4255 

+3.02 

-5.34 



15 

1.282 

1332 

“ 

4253 

+2.19 

-3.51 



16 

1.281 

1295 


4249 

-0.91 

+3.43 



17 

1.280 

1283 

“ 

4247 

+1.12 

-1.13 



18 

1.279 

1271 

“ 

4246 

-0.88 

+3.37 



Normal equations 



Solution 





1.33915 

-l.:37DOn 

+5.406 

= 0 

5 - -3V704 = 

“q 




-1.37005 +2.l4l9n 

-5 . 771 

= 0 

n *= 0.325 

dt Q 

+0 

'.'364 



[L 2 ] 

- 132. 

80 

tpvv] = 110.90 



Weight coefficients 



in 2 = 6.931 

m “ + 

2 V 

'63 


2.1608 

+1.3821 



6 ab = - 26 ° 48 

'03V70 

+ 

3 " 87 


+1.3821 

1.3509 



fc AB * - 75 ° 10 

+1 

vD 

o 

3V43 


Error ellipse 

Semiaxis 4771 

1, 48 

4. Adjustment of a Network 

For each station of the network, we need to calculate good approxi- 
mate coordinates designated as <j>, A, H. Using any reference ellipsoid 
with semiaxes a, b and denoting 



a 


we compute the Cartesian coordinates 


direction 37° 
127 


- 15 - 



x * (N + H) cos <j> cos A 

( y = (N + H) cos <j) sin A 

z = (N - e 2 N + H) sin <j> 


where 


I W = /l - e^ sin^<f> 


( 17 ) 


The fictitious observations 6 A _, t._ between two stations A and B 

AB AB 

give residuals 


/ 






Z B - Z A 


arc tan 


- 6 


/ 


AB 


<*B - X A )2 + (y B ‘ V 2 


(arc tan 


y B ' y A 


^ ' X A 


" t AB ) COS 6 AB 


(18) 


Differentiation of (18) shows that if x, y , z are corrected by dx , dy , dz, 
respectively, the new residuals will be 


(19) 

( 20 ) 

where 


v k = VdXg " dx A ) + B l (dy B " dy A ) + C l (dZ B " dz A ) + k AB 
I V h = A 2 <dX B - d V + B 2 <dy B ‘ dy A ) + h AB 
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A = sin 6 cos t 

1 r AB AB 


B = - 2. S i n 5 g i n t 

1 r AB AB 


C = + 2- cos 6 
1 r AB 


A 2 ' - r Bl " ‘ab 


B 2 ‘ + r cos ‘AB 


- / 


1 ‘ /(X B - X A >2 + <y B ' y A )2 + <Z B - Z A )2 


( 21 ) 


The weight matrix of (19) and (20) is given by (15). Therefore, we form 
equations 


(19) [paa] + (20) [pab] 
m /[paa] 


£21 _ iiakli . 

m V Lr [paa] 2 

\ 


( 22 ) 


which can be considered as uncorrelated observation equations with weights 
of 1. Then the normal equations can be formed and solved, in order to 
obtain corrections dx, dy, dz . 

In this adjustment, only the earth-fixed space directions have been 
used as observed quantities. Therefore, the network cannot be fixed with 
respect to any absolute zero point nor to the scale. If the corrections 
of all coordinates are considered to be unknown, the matrix of normal 
equations becomes singular. In fact, one point remains undetermined and, in 
addition, one coordinate of some other point. 

In the practice, the scale can be determined from the existing ground 
tri angulations between some stations. Perhaps, also, the observed distances 
of satellites from the ground stations can be used, with appropriate 
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weights. Similarly, the absolute zero point can be determined from existing 
astronomical observations of latitude and longitude, with appropriate 
gravimetric corrections for the deflection of the vertical. It is probable 
that the weights for all these additional observations are smaller than those 
for the space directions. Therefore, it is hardly feasible to include them 
in the main adjustment of the satellite network. 

In the following adjustments, the missing four coordinates are 
determined as follows: The observation equations (22) are taken in matrix 

form 


Ax 4* 1 * v 


(23) 


but the four coordinates cannot be determined. Therefore, we write 

x x + A 2 x 2 + 1 - v (24) 

and take four more equations 

B 1 X 1 + B 2 T x 2 = 0 (25) 

from conditions that the arithmetical mean of all ( or, at least, the best) 
approximate coordinates and of the side lengths will not be changed: 

[dx] = [dy] = fdz] = [xdx + ydy + zdz] « 0 (26) 


Because B^ is a 4 x 4 matrix, we can compute 



and (24) becomes 

(A 2 - Aj B^ 1 B 2 ) x 2 + 1 - v 

This system of equations is no longer singular. Therefore, denoting: 
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we can compute 



and the weight coefficients: 



Finally, we may point out that the Cartesian coordinates (17) are 
needed for the computation of residuals (18) only. The observation equa- 
tions (19) and (20) can be written using the original parameters <f>, A, 
and H in which purpose it is advisable to introduce temporary differentials 


dB - (M + H) d<j) 

dL - (N + H) cos <j> dA 

i 2 i 2 

1 - e 1 - e 

where M = ~ — N = r — a 

W W 


(27) 
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For the computation of the coefficients of these new observation equations 
the following matrix formula can be used: 



There are as yet no extensive observed nets which would provide data 
that could be used in numerical examples or in testing the methods 
suggested above. However, I am currently studying a schematic net which 
is a very good approximation of the global satellite network planned for 
international cooperation. In brief, I computed a net of 20 spherical 
triangles with 12 comer points and 30 sides, using the well-known proper- 
ties of a regular icosaeder. The middle points of each side were then 
taken as additional comer points. In this way, a network with 80 
triangles, 42 stations, and 120 sides was obtained. I hope that the 
experiences obtainable from the computation of such a regular system will 
give useful ideas to others for the planning of future computations of 
the networks actually observed. 
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